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摘要 : 球状 星团 和 星系 核 中 拥有 非常 密集 的 恒星 质量 黑洞 ， 当 两 个 黑洞 相互 靠近 时 ， 会 形成 
心率 极 高 的 轨道 并 发 射 引 力 波 。 该 类 信号 的 探测 将 有 助 于 黑洞 物理 的 进一步 研究 ， 力 

群 的 理解 。 然 而 ， 这 类 引力 波 尚 未 被 探测 到 ， 在 当前 的 研究 里 没有 被 广泛 考虑 。 因 此 ， 研 究 它 的 
可 探测 性 是 重要 的 。 根 据 数值 相对 论 模 拟 ， 最 近 的 研究 给 出 了 抛物 线 黑 洞 交会 的 准确 波形 。 利 用 
引力 波 搜索 工具 Minke， 在 总 质量 和 距离 两 个 参数 空间 里 制作 了 大 规模 的 数据 集 ， 并 定量 i 
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探测 范围 。 结 果 表 明 ， 黑 洞 交 会 产生 的 引力 波 很 可 能 是 当前 地 基 探 测 器 的 潜在 探测 对 象 ， 同 时 ， 
下 一 代 引 力 波 探测 器 明显 扩大 了 探测 范围 。 最 后 讨论 了 如 何 判 断 并 合 信号 为 黑洞 俘获 的 热点 问 
题 。 利 用 信号 在 两 种 采样 器 下 的 后 验 分 布 ， 提 出 了 一 种 区 分 方法 。 通 过 比较 不 同 的 后 验 ， 发 现 双 
黑洞 并 合 与 黑洞 俘获 信号 在 并 合 时 间 to 上 存在 显著 差异 ， 这 将 有 助 于 在 并 合 样本 中 区 分 这 两 种 


波源 。 

X 键 la: 引力 波 ; 恒星 级 黑洞 ， 双 黑洞 系统 
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1 引 言 


爱 因 斯 坦 在 1915 年 发 表 了 广义 相对 论 ， 并 预言 了 引力 波 ” 的 存在 ， 引 发 了 长 达 一 个 世 
纪 的 对 引力 波 信号 的 探索 。 这 项 工作 在 2015 年 9 月 获得 突破 ， 地 面 激光 干涉 仪 LIGO™ 组 
织 观测 到 由 一 个 双 黑 洞 并 合 产生 的 引力 波 信号 GW150914”， 为 天 体 和 字 宙 的 观测 提供 了 一 
种 新 的 手段 。 这 是 人 类 历史 上 首次 观测 到 引力 波 ， 但 是 LIGO 和 其 他 探测 器 的 使 命 还 远 未 
结束 ， 它 们 将 继续 探测 和 解释 探测 到 的 信号 。 引 力 波 探测 的 时 代 已 经 到 来 ， 随 之 而 来 的 是 解 
释 无 数 信号 ， 将 它们 与 噪声 区 分 开 ， 并 通过 分 析 深 入 了 解 产 生 这 些 引 力 波 的 天 体 物理 机 制 。 

在 过 去 的 几 年 里 ， 引 力 波 探测 器 探测 到 的 致密 双星 并 合 事件 急剧 增加 ，LIGO-Virgo” f& 
测 器 网 络 公布 的 引力 波 暂 现 源 星 表 1” 和 2" 包含 了 共计 55 个 引力 波 事件 。 在 最 新 公布 的 引 
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力 波 暂 现 源 星 表 3 中， 研究 小 组 又 新 确定 了 35 个 引力 波 事件 ， 其 中 有 32 个 最 有 可 能 是 
黑洞 并 合 ”"， 这 使 得 自 2015 年 来 探测 到 的 引力 波 事件 总 数 达到 90 个 ， 而 双 黑 洞 并 合 系统 已 
经 成 为 其 中 最 重要 的 引力 波源 。 其 形成 通道 引发 了 科学 家 的 巨大 兴趣 和 争论 ， 普 遍 认为 ， 恒 
星 级 双 黑 洞 并 合 的 形成 通道 主要 有 两 种 ， 一 种 是 孤立 的 大 质量 双星 演化 ， 另 一 种 则 是 两 
个 恒星 起 源 的 黑洞 在 动态 环境 中 形成 ””。 后 者 可 细 分 为 多 种 情形 ， 其 中 有 一 种 引力 波 

源 在 当前 的 搜索 中 并 未 被 广泛 地 考虑 ， 即 黑洞 交会 (black hole encounter) ™ 
球状 星团 和 星系 中 心 的 空间 区 域内 存在 非常 密集 的 致密 星 ， 为 黑洞 之 间 的 相互 作用 提 
供 了 一 个 活跃 的 环境 。 这 里 存在 着 多 种 引力 作用 过 程 : 其 中 一 些 黑洞 双星 形成 束缚 轨道 ， 发 
射出 具有 旋 进 、 并 合 与 铃 宕 三 个 阶段 的 周期 性 引力 波 ， 在 旋 进 过 程 中 引力 波 带 走 能 量 使 轨道 
圆 化 ， 从 而 演化 为 圆 轨道 或 低 偏心 率 椭 圆 轨 道 的 双星 系统 ， 最 终 产生 并 合 事件 ， 这 类 系统 被 
黑洞 并 合 (binary black hole merger); 而 其 他 的 黑洞 双星 则 会 形成 非 束 缚 轨道 ， 导 致 
能 的 散射 或 并 合 事件 。 非 束缚 轨道 的 形成 条 件 是 : 

bob 之 0 ， (1) 
> 其 中 ，E 是 双星 系统 的 轨道 能 量 。 我 们 将 非 束缚 轨道 的 黑洞 双星 称 为 黑洞 交会 ， 它 们 发 
<+ 射 强烈 的 爆发 型 引力 波 。 与 双 黑 洞 并 合 相 比 ， 其 波形 没有 旋 近 阶段 ， 且 进入 LIGO 探测 器 有 
效 波段 时 具有 非常 显著 的 偏心 率 (90%， 偏 心率 < > 0.9)， 在 物理 情景 上 可 以 与 前 者 区 分 开 
来 。 在 轨迹 偏转 角 较 小 的 情况 下 ， 这 个 过 程 可 以 被 认为 类 似 于 电磁 辐射 中 的 轨 致 加 射 宇 汪 ， 
但 是 当 考 虑 到 较 大 的 偏转 角 和 自 旋 时 ， 引 力 辐 射 变 得 更 加 复杂 。 徙 致 辐射 情况 下 的 黑洞 交会 
波形 已 经 提出 了 很 好 的 近似 模型 来 解释 ， 它 适用 于 具有 任意 轨道 偏转 角 的 低速 情况 ”和 黑 
洞 的 正面 碰撞 中。 根据 轨道 能 量 的 条 件 ， 数 值 相对 论 建 模 得 到 的 黑洞 交会 波形 又 可 由 轨 
道 类 型 分 为 双 曲 线 型 (Bw > 0) 和 抛物 线 型 (Eob = 0)。 最 近 的 数值 模拟 进展 已 经 允许 为 无 
自 旋 黑洞 双星 的 双 曲 线 交 会 产生 3.5 fr PNE 的 波形 ”， 为 无 自 旋 黑 洞 双星 的 抛物 线 交 会 产 
生 准 确 的 波形 。 进一步 地 ， 当 两 个 黑洞 相互 靠近 时 ， 如 果 引 力 波 辐 射 带 走 的 能 量 超过 轨 
道 能 量 ， 系 统 仍 将 形成 束缚 轨道 并 产生 并 合 事件 ， 部 分 的 黑洞 交会 称 之 为 黑洞 俘获 (black 

hole capture); 反之 ， 双 星 会 发 生 散射 ， 经 过 引力 的 相互 作用 ， 先 相互 靠近 而 后 永远 分 开 。 
关于 这 类 波源 ， 除了 前 文 所 到 的 波形 模拟 ， 现 有 的 研究 还 关注 其 能 谱 的 计算 ”和 探 
测 率 的 估计 。 黑 洞 交 交会 事件 的 探测 率 可 以 限制 星系 核 中 的 恒星 初始 质量 函数 和 中 心 黑洞 质 
量 ~”， 有 几 项 研究 已 对 先进 LIGO Al LISA 中 此 类 事件 的 单位 距离 探测 率 作出 估计 ， 

其 结果 都 为 每 年 约 1 Gpe-s 的 数量 级 但 这 类 结果 仍 缺 乏 对 引力 波 信号 探测 范围 的 估计 。 
本 文 使 用 抛物 线 黑洞 交会 波形 ， 分 别 对 并 合 、 非 并 合 两 种 情况 的 波形 制作 模拟 数据 集 ， 
并 对 多 个 探测 器 下 的 探测 范围 展开 研究 。 最 新 的 地 面 和 空间 引力 波 探测 器 可 能 探测 到 的 黑洞 
交会 引力 波 是 研究 的 主要 动机 。 由 于 黑洞 交会 缺乏 良好 的 模板 和 搜索 算法 ， 在 先进 LIGO, 
先进 Virgo 和 KAGRA 组 成 的 探测 器 网 络 的 前 三 轮 观测 周期 里 并 未 有 类 似 的 事件 得 到 探测 
和 确认 。 未 来 投入 使 用 的 LISA、DECIGO 和 爱 因 斯 坦 望 远 镜 (Einstein Telescope, ET)， 将 


?post-Newtonian (PN): 对 于 广义 相对 论 中 引力 场 的 微 扰 解 ， 在 2o« 1 时 ， 将 其 用 级 数 展开 ，z 次 方 数 的 1/2 叫做 
PN 的 阶 数 。 
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加 入 探测 器 网 络 并 贡献 信 噪 比 更 高 的 观测 数据 ， 可 以 把 引力 波 的 探测 频段 拓展 到 更 低频 、 更 
宽 的 范围 ， 因 此 研究 这 类 事件 的 可 探测 性 很 有 意义 。 

关于 该 类 引力 波源 中 的 并 合 波形 ， 另 一 个 值得 研究 的 问题 是 如 何 将 并 合 的 波形 与 探测 
到 的 高 质量 双 黑 洞 并 合 信号 正确 区 分 ”。 如 前 所 述 ， 典 型 的 双 黑 洞 并 合 引 力 波 具 有 旋 进 、 
并 合 与 铃 宕 三 个 阶段 ， 其 波形 形状 与 爆发 式 的 黑洞 俘获 引力 波 非 常 不 同 。 然 而 ， 在 总 质量 
很 高 的 情形 下 ， 前 者 波形 中 的 “ 旋 进 ”阶段 低 于 LIGO 探测 器 的 截止 频率 ， 不 能 反映 到 记 
录 数 据 里 ， 因 此 这 两 类 事件 会 具有 相似 的 波形 ， 在 实际 的 工作 中 可 能 被 误 分 类 。 典 型 的 例 
子 如 GW190521™ 四 ， 它 由 一 个 双 黑洞 系统 并 合 发 出 ， 两 个 黑洞 的 质量 分 别 为 85 Mo 和 66 
Mo， 并 合 得 到 了 一 个 150 Mo 的 中 等 质量 黑洞 多 个 研究 ”表明 ， 它 的 轨道 可 能 具有 很 训 
的 偏心 率 ， 因 此 引发 了 对 其 形成 机 制 的 关注 。 最 近 的 研究 已 证 明 ， 将 黑洞 俘获 用 主流 的 圆 轨 
道 双 黑 洞 并 合 模板 IMRPhenomPv2”” 进行 分 析 ， 仍 可 得 到 类 似 高 质量 双 黑 洞 并 合 的 贝 叶 
斯 推理 结果 ， 这 对 引力 波 事件 的 正确 识别 和 分 类 提出 了 挑战 。 
本 文 介绍 了 黑洞 交会 ， 一 种 潜在 的 引力 波源 ， 预 计 将 在 新 一 代 地 基 引 力 波 探测 器 和 空间 
引力 波 探测 器 中 具有 可 观 的 探测 率 。 本 文 还 利用 数值 相对 论 模 拟 的 波形 ， 在 多 个 参数 空间 里 
创建 了 该 类 波源 的 大 规模 引力 波 数据 集 ， 然 后 通过 信 噪 比 的 计算 工作 揭示 了 探测 范围 和 总 
质量 以 及 光度 距离 的 关系 ， 并 探讨 了 其 可 探测 性 。 同 时 ， 我 们 还 将 下 一 代 探 测 器 爱 因 斯 坦 望 
远 镜 的 探测 效果 与 先进 LIGO 进行 比较 。 最 后 ， 针 对 黑洞 交会 和 高 质量 双 黑 洞 并 合 这 两 种 
引力 波源 在 实际 数据 分 析 工 作 中 可 能 被 混淆 的 问题 ， 我 们 提出 了 一 种 方法 以 区 分 。 对 黑洞 交 
会 信号 的 两 种 后 验 概率 分 布 ， 我 们 计算 了 它们 在 几 个 主要 参数 上 的 差别 ， 并 用 JS 散 度 衡量 
其 大 小 ， 然 后 将 它们 与 双 黑 洞 并 合 的 基准 值 进行 了 比较 。 

本 文 的 结构 如 下 ， 第 2 章 介 绍 了 生成 大 规模 黑洞 交会 引力 波 波形 的 方法 ， 利 用 最 新 的 
10 个 数值 相对 论 波形 ， 包 括 5 个 并 合 波形 和 5 个 非 并 合 波形 ， 来 制作 模拟 数据 集 ， 进 而 
算 它 们 的 信 噪 比 和 可 探测 性 。 第 3 章 提出 了 一 种 区 分 双 黑 洞 并 合 和 黑洞 交会 的 方法 ， 通 
比较 两 种 采样 器 输出 的 后 验 分 布 ， 对 4 种 抛物 线 黑洞 俘获 引力 波 波形 进行 分 析 。 第 4 章 是 
总 结 和 展望 ， 指 出 了 现 有 工作 中 可 以 提高 的 部 分 ， 并 对 未 来 的 工作 方向 进行 了 讨论 。 


im 
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= 2 可 探测 性 估计 


黑洞 交会 是 一 种 尚未 被 探测 到 的 引力 波源 。 在 本 章 中 ， 我 们 使 用 抛物 线 黑洞 交会 波 
制作 模拟 数据 ， 以 研究 它们 在 当前 和 未 来 探测 器 中 的 可 探测 性 。 
2.1 ”模拟 数据 注入 

引力 波 的 探测 器 应 变 h(t) RA hy 和 hh 两 个 分 


NS 


fan 


A(t) = F (t)h4 (t) + Fx(0)h«(t) ， (2) 

1 
F= 5( + sin? 6) cos 2a cos 24) — sindsin2asin2y , (3) 
Fg = jü + sin? 6) cos 2a sin 24) — sind sin2acos2y , (4) 
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Hp, Fy PE. 是 分 别 对 应 于 引力 波 两 个 极 化 分 量 的 天 线 响应 ，6 是 引力 波源 相对 于 探测 
器 平面 的 高 度 角 ，a 是 相对 于 其 中 一 个 辟 的 方位 角 ， 而 是 引力 波 的 极 化 角 。 由 此 可 见 ， 
源 的 天 空 方位 角 将 对 信号 的 波形 及 幅度 产生 影响 ， 而 天 线 响 应 还 会 降低 探测 器 对 天 空 某 些 
区 域 的 灵敏 度 。 为 了 计算 最 佳 信 噪 比 ， 我 们 在 随后 的 数据 模拟 中 选取 合适 的 天 空 方位 角 ， 以 
保证 探测 器 的 探测 效果 处 于 最 佳 水 平 。 

我 们 使 用 Minke 软件 包 来 产生 模拟 引力 波 信号 。Minke 是 一 个 基于 Python 语言 的 引力 
波 搜 索 工具 ， 可 以 根据 指定 的 数值 相对 论 波形 和 注入 参数 分 布 制作 大 规模 的 数据 集 。 使 用 
Minke 进行 产生 模拟 数据 的 过 程 如 下 。 

(1) 指定 数值 相对 论 波形 。 本 研究 使 用 抛物 线 黑洞 交会 波形 ， 包 括 并 合 类 和 非 并 合 类 ， 
每 类 又 根据 质量 比分 为 1, 2, 4, 8, 16 五 种 ， 总 计 有 10 种 波形 。 


4 


u 


(2) 设置 注入 参数 的 分 布 ， 包 括 黑洞 双星 的 总 质量 、 源 光度 距离 、 天 空 方位 角 、 源 的 极 
化 角 。 
(3) 调用 LALSimulation 为 每 一 组 注入 参数 产生 模拟 波形 
T (4) 指定 探测 器 ， 根 据 其 天 线 方向 图 计算 探测 器 对 每 一 个 模拟 波形 的 天 线 响应 ， 得 到 模 


拟 信号 集 。 这 里 分 别 使 用 了 先进 LIGO， 爱 因 斯 坦 望远镜 和 DECIGO。 

(5) 最 终 ， 为 每 一 个 引力 波 信号 注入 探测 器 噪声 ， 获 得 模拟 数据 集 。 

我 们 设置 注入 的 双星 总 质量 为 具有 50 个 点 的 均匀 分 布 ， 范 围 为 5 Mo ~ 150 Mo. $% 
洞 根据 质量 可 以 分 为 恒星 级 黑洞 、 中 等 质量 黑洞 和 超大 质量 黑洞 ， 其 中 ， 已 知 的 中 等 质量 黑 
洞 极其 少 ， 黑 洞 质量 分 布 在 恒星 级 黑洞 与 超大 质量 黑洞 之 间 存 在 间隙 。 由 于 探测 器 频率 的 
限制 ， 目 前 观测 到 的 并 合 事件 都 是 来 自 两 个 恒星 级 黑洞 的 双星 系统 。 因 此 ， 我 们 采取 的 黑洞 
总 质量 上 限 为 150 Me， 它 来 源 于 迄今 为 止 观测 到 的 总 质量 最 大 的 并 合 事件 GW190521。 在 
这 个 质量 范围 里 ， 单 个 黑洞 为 宇宙 中 大 量 存 在 的 恒星 级 质量 黑洞 ， 且 波源 有 间接 的 光学 观 
测 可 相互 验证 。 我 们 设置 注入 的 光度 距离 的 最 小 值 为 1 Mpc， 根 据 每 个 波形 来 确定 不 同 的 
最 大 值 ， 也 形成 50 个 注入 值 的 均匀 分 布 。 信 和 号 的 起 始 GPS 时 间 设 为 太 = 1126 259 642.0 s. 
我 们 先 对 先进 LIGO 注入 信号 ， 由 前 面 的 质量 和 距离 分 布 ， 每 一 个 波形 得 到 了 2500 个 信号 


(C 样本 。 图 四 展示 了 其 中 一 例 信号 的 频谱 特征 应 变 与 探测 器 灵敏 度 的 比较 ， 该 信号 由 “black 
= hole encounter” 表 示 ， 它 来 自 质量 比 g = 1 的 黑洞 俘获 波形 ， 总 质量 为 150 Me， 光度 距离 
"T 985000 Mpc. FAD RE T 10 种 黑洞 交会 波形 的 各 一 例 信号 ， 它 们 的 区 别 在 于 双星 质量 比 g 
~ 的 不 同和 双星 最 终 是 否 并 合 。10 种 波形 都 拥有 150 Mo 的 总 质量 ， 和 50 Mpe 的 光度 距离 。 
相 比 最 终 没有 并 合 的 黑洞 交会 ， 黑 洞 俘获 在 波形 上 具有 并 合 与 铃 宕 两 个 阶段 的 独特 形状 。 
2.2 _ 信 了 噪 比 限制 探测 范围 
引力 波 信号 的 最 佳 信 品 比 pop, 可 表示 为 : 
fmax |p2 2 
[Eg] : 


其 中 ，h(f) 是 引力 波 信 号 的 频 域 形式 ，S;(f) 是 以 Hz! 为 单位 的 单 边 噪声 谱 密 度 ，fiin < 
f S fmax 对 应 于 仪器 的 工作 频段 ， 例 如 先进 LIGO 有 (fmin, fmax) © (10,104) Hz。 可 以 看 
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应 变 


应 变 


“~~ ALIGO 
^ Tnitial LIGO 
“~~ Decigo 

— ET 


== black hole encounter 


特征 应 变 /Hz0s 


10-26L_ 
10? 107 10° 10! 10? 10° 104 
频率 /Hz 


iE: black hole encounter， 质 量 比 g = 1 的 黑洞 交会 波形 ， 总 质量 为 150 Mo ， 光 度 距 离 为 5000 Mpc. 


图 1 黑洞 交会 的 特征 应 变 与 多 种 引力 波 探 测 器 灵敏 度 的 对 比 
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Hoa) 黑洞 俘获 波形 ，b) 非 并 合 的 黑洞 交会 波形 。 质 量 比 用 q 表示 。 可 以 看 出 ， 两 种 波形 都 没有 旋 近 阶段 。 黑 
洞 俘获 波形 拥有 明显 的 并 合 、 铃 宕 阶段 ， 该 阶段 的 时 间 尺 度 为 0.1 s; 其 中 质量 比 为 2 的 波形 在 并 合 前 具有 额外 
的 振荡 。 黑 洞 交会 波形 则 只 有 短暂 的 振荡 过 程 ， 黑 洞 双星 由 引力 作用 相互 靠近 而 后 永远 分 开 。 


图 2 不同 质量 比 的 黑洞 交会 波形 和 黑洞 俘获 波形 
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出 ， 最 佳 信 噪 比 考 虑 了 噪声 功率 谱 密度 的 影响 ， 因 此 噪声 功率 谱 密 度 高 的 频段 对 信 噪 比 的 贡 
献 小 于 低 的 频段 。 

对 于 每 种 波形 ， 我 们 使 用 公式 (5) 为 每 个 信号 计算 信 品 比 ， 然 后 将 它们 绘制 成 图 ， 如 图 
加 和 图 团 所 示 。 


201 — EE 
2000 4000 6000 8000 10000 1000 3000 5000 7000 
SR ES /Mpc 距离 /Mpc 
a) b) 

PA 140 MHH / 
> | 
= 100; 
ag: 
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距离 /Mpc 距离 /Mpc 
c) d) 
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= 100} z 

hr 

"E 


20 (AA 
X 300 400 600 800 1000 
SR ES /Mpc 
e) 

YE: a)q=1; b)q=2; c)g=4; d)q=8; e) g = 二 16。 信 品 比 为 8 的 曲线 代表 了 具有 该 信 噪 比 的 所 有 信 
号 样本 。 位 于 该 曲线 右 下 方 的 信号 样本 ， 等 值 线 梯度 下 降 ， 其 信 噪 比 均 低 于 8， 因 此 不 能 被 先进 LIGO 探测 到 ; 
位 于 该 曲线 左上 方 的 信号 样本 由 于 信 品 比 高 于 8 而 可 以 被 探测 。 由 于 总 质量 的 上 限 限制 ，5 个 波形 都 在 纵 坐 标 
(双星 总 质量 ) 最 大 时 取得 最 大 探测 距离 ， 且 黑洞 双星 的 质量 比 g 越 接近 1， 其 产生 的 引力 波 可 探测 范围 越 大 。 


3 FË LIGO 中 抛物 线 黑洞 俘获 波形 的 信 噪 比 等 值 线 图 


TN 


E] BL Jos SAP RTE RU fS ERES. “EG BER PLNS JAY ER IC. fer EE BÉ i EE BB 
(AA 8, fa SWAB BE LACER SI, JF BOAR. 3 清楚 地 表明 ， 来 自 这 
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e) 


i: a) q=1; b) q = 2; c) q = 4; d) q = 8; e) q = 16. HR 3, 不同 的 是 ， 在 恒星 质量 黑洞 的 范围 内 ， 
信 噪 比 等 值 线 出 现 拐点 ， 此 拐点 横 坐 标 便 对 应 波形 的 最 大 探测 距离 ， 且 黑洞 双星 的 质量 比 q 越 接近 1， 取 得 拐点 
的 纵 坐 标 (双星 总 质量 ) 越 小 。 


4 ”先进 LIGO 中 非 并 合 抛物 线 黑洞 交会 波形 的 信 品 比 等 值 线 图 


样 一 个 系统 的 引力 波 ， 在 双星 总 质量 最 大 为 150 Mo 的 假设 下 ， 最 远 可 以 在 6500 Mpc ( 质 
量 比 为 1 的 黑洞 双星 ) 的 距离 探测 到 ， 考 虑 到 联合 观测 对 灵敏 度 的 提升 ， 这 个 数值 可 能 会 
更 高 。 同 时 ， 质 量 比 为 2, 4, 8, 16 的 双星 分 别 具 有 5500, 1800, 1250, 580 Mpc 的 最 大 探测 
距离 ， 质 量 比 越 接近 于 1 的 系统 探测 范围 越 大 。 图 团 展 示 了 非 并 合 黑 洞 交会 波形 的 信 噪 比 。 
最 大 探测 距离 在 20 Mo ~ 60 Mo 的 总 质量 范围 内 取得 ， 分 别 为 57, 118, 115, 275, 160 Mpc, 
对 应 质量 比 为 1, 2, 4, 8, 16 的 黑洞 双星 。 因 此 ， 黑 洞 俘获 比 非 并 合 的 黑洞 交会 在 先进 LIGO 
中 的 探测 范围 更 远 ， 也 更 有 和 希望 被 发 现 。 
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在 考虑 下 一 代 探 测 器 对 黑洞 交会 引力 波 的 探测 能 力 时 ， 我 们 使 用 爱 因 斯 坦 望 远 镜 和 
DECIGO 探测 器 执行 了 与 先进 LIGO 相同 的 信号 注入 与 信 噪 比分 析 。 结 果 表 明 ， 对 于 
黑洞 俘获 ， 爱 因 斯 坦 望 远 镜 在 双星 质量 比 g = 1 的 情况 下 获得 最 大 探测 距离 ， 可 达到 
5x10* Mpc; 对 于 非 并 合 的 黑洞 交会 ， 则 可 以 在 最 远 2500 Mpc 的 距离 上 被 探测 到 ， 此 时 对 
应 gq = 8 的 波形 。 另 一 方面 ，DECIGO 对 黑洞 俘获 的 探测 效率 一 般 ， 最 大 可 达 3000 Mpc, 
氏 于 目前 正在 运行 的 先进 LIGO 的 探测 能 力 ; 但 其 对 非 并 合 的 黑洞 交会 拥有 最 大 的 探测 范 
围 ， 最 远 可 达 5000 Mpc， 此 时 波形 质量 比 有 g=1 ( 


30 000 50 000 


5 000 


SR ES /Mpc 距离 /Mpc 
c) d) 


YE: a) q = 1 的 黑洞 俘获 ， 爱 因 斯 坦 望 远 镜 ; b) g = 1 的 黑洞 俘获 ，DECIGO; c) — 8 的 非 并 合 黑洞 交 
会 ， 爱 因 斯 坦 望远镜 ; d) 9 = 1 的 非 并 合 黑洞 交会 ，DECIGO。 同 图 3， 不 同 的 是 ， 本 文 以 爱 因 斯 坦 望 远 镜 和 
DECIGO 为 例 ， 分 别 选 取 了 黑洞 俘获 和 非 并 合 黑洞 交会 里 拥有 最 大 探测 范围 的 等 值 线 图 作为 展示 。 


图 5 下 一 代 探测 器 中 抛物 线 黑洞 交会 波形 的 信 噪 比 等 值 线 图 


3 利用 JS 散 度 区 分 黑洞 俘获 与 双 黑 洞 并 合 


目前 的 参数 估计 工作 主要 采用 双 黑 洞 并 合 模 板 IMRPhenomPv2 执行 贝 叶 斯 推理 。 根 据 
贝 叶 斯 理论 ， 后 验 概率 分 布 可 表示 为 : 

p(xly, H) x pyle, H)p(zx. H) ， (6) 
其 中 ，z 是 模型 参数 ，%y 是 观测 数据 ， 互 RAY, p(vx|y,H) JESUS. p(y|v,H) 是 似 然 函数 ， 
而 pls, H) 是 这 些 参数 的 先 验 。 在 公式 右边 我 们 省 略 了 一 个 比例 系数 ， 贝 叶 斯 证 据 p(y) 的 
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倒数 ， 这 么 做 是 因为 p(y) 是 常数 且 我 们 只 对 后 验 分 布 的 形状 感 兴趣 。 这 类 分 析 对 双 黑 洞 并 
合 与 黑洞 俘获 的 推理 存在 简 并 性 ”。 我 们 从 最 近 的 研究 ” 中 获得 黑洞 俘获 和 双 黑 洞 并 合 
事件 主要 参数 的 后 验 概 率 分 布 ， 它 们 在 双 黑 洞 并 合 的 零 自 旋 IMRPhenomPv2 近似 下 使 用 
dynesty fll Vitamin” 两 种 采样 器 分 别 进行 了 参数 估计 。dynesty 采用 了 嵌 套 采样 的 方法 ， 
而 Vitamin 基于 深度 学 习 技术 给 出 双 黑 洞 并 合 的 后 验 估计 。 
我 们 提出 一 种 基于 Js 散 度 思 的 区 分 方法 ， 用 来 在 实际 探测 中 判断 类 似 双 黑洞 并 合 的 信 
号 是 否 为 黑洞 俘获 。 这 一 方法 的 基础 思想 是 ， 一 个 模型 对 符合 其 物理 机 制 的 数据 具有 更 好 的 
推理 结果 。 在 双 黑 洞 并 合 模型 下 ， 该 类 信号 本 身 可 以 得 到 非常 精确 的 推断 ， 而 黑洞 交会 信号 
则 得 到 了 具有 偏见 的 推断 。 在 后 一 种 情况 下 ， 推 理 的 结果 仅仅 是 数学 的 ， 对 应 到 每 个 参数 上 
则 不 具有 实际 意义 。 因 此 ， 当 我 们 比较 两 种 采样 器 的 偏见 推断 时 ， 每 个 参数 的 后 验 分 布 可 能 
了 很 大 的 差异 。 而 对 于 符合 物理 机 制 的 信号 ， 不 同 采样 器 的 推断 将 趋 近 于 一 致 。 我 们 用 JS 
散 度 来 量化 这 种 差异 ， 则 对 于 一 个 具有 类 似 双 黑洞 并 合 后 验 的 未 知 信号 ，JS 散 度 的 大 小 便 
可 以 表明 信和 号 是 否 来 源 于 该 模型 的 物理 机 制 。 这 种 方法 提供 了 一 种 判 据 ， 可 以 从 双 黑 洞 并 合 
- 事件 中 找 出 如 黑洞 俘获 之 类 的 未 建 模 或 不 良 建 模 的 瞬 态 信号 。 
JS 散 度 是 衡量 两 个 分 布 p(z) 与 g(z) 之 间 差异 度 的 量 ， 其 定义 为 : 


Dis(p| d) =3[Drr(p|s) + Dr(a l )] ， (7) 
Diap| à) = f ooog (22) ar , (8) 


KH, s=1/2(p+q), Der 是 KL 散 度 。 
JS 散 度 在 [0, 1] 取 值 , 其 值 越 大 说 明 两 个 分 布 的 差异 越 大 。 接 下 来 ， 我 们 考虑 这 种 分 

析 方 法 选取 的 参数 。 无 自 旋 的 圆 轨道 双 黑 洞 并 合 模型 具有 9 个 参数 ， 分 别 是 双 黑 洞 的 质量 
(mi, m) EREB (dL)、 天 空 方位 角 (o, 6)、 双 星 的 轨道 倾角 (Omn) IRRA (v). 
并 合 时 间 (tm). HAHM (bo)。 选 取 参 考 时 间 tret = 1126 259642.5 s， 则 相对 并 合 时 间 为 
to = tm — tep RIDER mi, mo, dr 和 to 分 别 进行 散 度 分 析 。 对 于 dynesty 和 VItamin X 
样 器 输出 的 后 验 分 布 ， 我 们 先 对 两 种 后 验 中 的 样本 数 较 大 的 一 个 进行 随机 抽样 ， 获 得 的 新 样 
本 和 样本 数 较 小 的 后 验 一 致 ， 然 后 对 m, ma, di, Fl to 4 个 参数 的 分 布 分 别 进行 JS 散 度 的 
T 计算 和 比较 。 我 们 选取 了 4 种 质量 比 的 抛物 线 黑洞 俘获 波形 (q = 1, 4, 8, 16) 和 一 例 零 自 旋 
的 双 黑 洞 并 合 信号 执行 这 种 分 析 。 篇 幅 所 限 ， 这 里 的 黑洞 俘获 仅 展 示 了 质量 比 为 1 的 波形 
的 后 验 分 布 ， 如 图 日 所 示 ; 同时 还 展示 了 双 黑 洞 并 合 波形 作为 对 比 标准 ， 如 图 四 所 示 。 所 

的 JS 散 度 值 记录 在 表 加 中 。 
从 图 加 和 图 四 可 以 看 出 ， 与 dynesty 相 比 ，VItamin 对 两 种 波形 里 m, 和 ma 的 采样 都 
具有 更 尖锐 的 峰 。 此 外 ，dynesty 和 Vitamin 对 双 黑 洞 并 合 的 并 合 时 间 to 具有 非常 相近 的 
估计 ， 对 黑洞 俘获 则 不 具有 这 种 特点 ，VItamin 的 to 后 验 峰 值 比 真 值 和 dynesty 的 后 验 峰 
值 晚 了 约 0.03 s。 这 点 从 表 四 反映 得 更 为 清楚 ， 两 者 的 JS 散 度 分 别 为 0.05 和 0.22。 除 了 质 
量 比 为 16 的 波形 ， 其 他 黑洞 俘获 也 获得 了 类 似 的 结果 ，VItamin 采 样 器 总 是 倾向 于 把 如 估计 
为 具有 双 峰 ， 且 主峰 远离 真 值 的 后 验 分 布 ， 这 导致 如 具有 很 高 的 JS 散 度 。 质 量 比 为 16 的 
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m,/Mo m,/ Mo 


注 : Vitamin Æ m, fll ma 上 上 有 具有 更 陡峭 的 峰 ; 
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两 种 采样 器 在 to 上 具有 位 置 相同 的 双 峰 ， 但 主峰 位 置 不 一 致 。 


6 ”抛物线 黑 洞 俘获 的 后 验 分 布 比较 ， 使 用 dynesty 和 Vitamin 两 种 采样 器 


表 1 dynesty 和 VItamin 采样 器 后 验 的 JS RE 


mai m2 dL to 
q=1 0.10 0.16 0.03 0.22 
q=4 0.32 0.05 0.19 0.36 
q=8 0.05 0.28 0.14 0.34 
q=16 0.27 0.16 0.08 0.06 
双 黑 洞 并 合 0.19 0.09 0.07 0.05 
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图 7 


YE: VItamin Æ mı M ma 上 具有 更 陡峭 的 峰 ; 


双 黑 洞 并 合 的 后 验 分 布 比较 ， 使 用 dynesty 和 Vitamin 两 种 采样 器 


两 种 采样 器 


d 


1600 2400 0.216 0.219 0.222 
d,/Mpc t,/s 


在 to 上 具有 单 峰 。 


波形 在 如 上 的 散 度 较 小 ， 这 是 因为 它 的 两 个 后 验 主峰 几乎 重合 。 在 其 他 几 个 参数 上 ， 黑 洞 


俘获 的 JS 散 度 则 没有 明显 : 


也 与 双 黑 洞 并 合 区 分 开 


散 度 分 析 为 区 分 双 黑 洞 


并 合 与 黑洞 俘获 提供 


来 。 
了 一 种 判 据 : 当 一 个 类 似 双 黑洞 并 合 的 未 


知 信号 被 探测 到 时 ， 用 dynesty 和 Vitamin 分 别 获得 并 合 时 间 如 的 后 验 分 布 ， 并 计算 它们 
之 间 的 JS 散 度 ， 若 其 值 远大 于 双 黑 洞 并 合 具 有 的 值 0.05， 则 表明 其 很 有 可 能 是 黑洞 俘获 


信号。 


4 总 结 与 展望 


本 文 讨论 了 抛物 线 黑 洞 交会 引力 波 的 探测 。 本 研究 的 主要 目标 有 两 个 方面 : 
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(1) 在 总 质量 和 距离 两 个 参数 空间 上 产生 大 批量 的 引力 波 信号 样本 ， 用 以 评估 先进 
LIGO 和 下 一 代 引 力 波 探测 器 对 多 个 波形 的 探测 能 
(2) 提出 一 种 区 分 方法 ， 通 过 比较 信和 号 的 两 种 后 验 分 布 来 解决 黑洞 俘获 和 高 质量 双 黑 洞 
并 合 在 实际 工作 中 的 后 验 结果 简 并 问题 。 
第 一 项 研究 的 结果 表明 ， 在 双星 最 大 总 质量 为 150 Mo 的 前 提 下 ， 目 前 的 先进 LIGO 可 
以 对 最 远 6500 Mpe 的 黑洞 交会 产生 可 探测 的 信 噪 比 。 值 得 注意 的 是 ， 即 使 是 不 同 波形 中 
最 小 的 探测 范围 ， 其 值 也 达到 了 57 Mpc， 超 过 本 超星 系 团 的 直径 ， 使 得 这 类 事件 成 为 非常 
希望 的 引力 波源 候选 者 。 先 进 LIGO 对 并 合 与 非 并 合 的 两 类 波形 取得 最 大 探测 范围 的 条 
件 也 不 同 。 黑 洞 俘获 的 总 质量 越 大 ， 其 探测 距离 越 大 ， 而 非 并 合 黑洞 交会 则 在 总 质量 为 (20 
~ 60) Mo 时 有 最 大 探测 距离 。 我 们 还 将 波形 注入 到 爱 因 斯 坦 望 远 镜 和 DECIGO 中 ， 比 较 
了 未 来 一 代 探 测 器 对 该 类 事件 探测 能 力 的 提升 。 爱 因 斯 坦 望 远 镜 将 大 大 提高 黑洞 俘获 的 最 
远 探测 距离 ， 最 远 可 达 5 x 104 Mpc, DECIGO 则 拓展 了 非 并 合 黑洞 交会 的 探测 范围 ， 最 远 
可 达 5000 Mpc。 在 第 二 章 ， 我 们 为 了 解决 参数 估计 工作 中 黑洞 交会 的 误 分 类 问题 ， 提 出 了 
= 一 种 基于 JS 散 度 的 区 分 方法 。 我 们 比较 了 波形 在 两 种 采样 器 Vitamin 和 dynesty 下 的 后 验 
Z 分 布 ， 这 些 后 验 分 布 数据 来 自 先前 的 一 项 研究 ” ， 包 括 4 个 黑洞 俘获 波形 和 1 个 双 黑洞 并 
5 合 波形 。 分 析 结 果 显 示 ， 黑 洞 俘获 在 并 合 时 间 to 这 一 项 拥有 的 JS 散 度 远 高 于 圆 轨 道 双 黑 洞 
并 合 ， 因 此 它 将 有 助 于 在 探测 到 的 并 合 样本 中 区 分 这 两 种 引力 波源 。 
根据 一 项 最 新 的 估计 ， 黑 洞 交会 的 事件 率 约 为 每 年 0.9 Gpc-3”。 本 文 计算 的 探测 范围 
可 以 表明 ， 黑 洞 交 会 的 事件 率 可 能 是 相当 可 观 的 ， 搜 索 它 们 是 对 已 有 观测 资源 的 充分 有 效 
利用 ， 并 且 将 在 新 一 代 探 测 器 中 持续 进行 。 为 了 在 引力 波 数据 中 提取 出 准确 的 信号 ， 我 们 需 
要 更 加 精确 的 相对 论 模拟 波形 ， 这 项 工作 可 能 会 促进 这 些 搜索 。 此 外 ， 对 于 第 二 项 研究 ， 我 
们 将 在 本 文 的 基础 上 进行 更 多 的 测试 和 验证 。 首 先 ， 受 到 计算 资源 的 限制 ， 目 前 所 用 的 事件 
略 少 ， 未 来 的 工作 将 考虑 更 多 模拟 波形 在 并 合 时 间 如 上 的 表现 。 例 如 ， 将 双 黑 洞 并 合 波形 
拓展 到 不 同 质量 比 ， 以 及 将 黑洞 俘获 波形 注入 不 同 的 总 质量 。 其 次 ， 本 文 使 用 的 后 验 分 布 来 
自 零 自 旋 的 双 黑 洞 并 合 模型 IMRPhenomPv2， 在 自 旋 、 轨 道 偏心 率 和 进 动 三 个 方面 可 以 进 
行 改进 。 通 过 动态 俘获 形成 的 黑洞 双星 很 可 能 具有 各 向 同性 的 自 旋 方 向 ”， 这 个 特征 是 孤 
立 大 质量 双星 演化 形成 的 双 黑 洞 系统 所 不 具备 的 ; 因此 ， 自 旋 的 分 析 对 于 这 类 事件 的 区 分 
有 具有 重要 意义 。 此 外 ， 该 模型 利用 单 自 旋 的 后 牛顿 旋转 来 描述 进 动 效应 ， 而 最 新 的 模型 允许 
= 双 自 旋 旋 转 ， 该 模型 还 只 限于 圆 轨道 的 双 黑 洞 系统 ， 因 此 可 以 考虑 在 下 一 步 引 入 偏心 轨道 
的 SEOBNRE 模型 ， 虽 然 它 的 建 模 也 无 法 达到 这 类 波源 的 超 高 轨道 偏心 率 ， 但 预计 有 比 贺 
轨道 模型 更 好 的 拟 合 效果 。 第 三 ， 主 流 的 参数 估计 工作 普遍 使 用 了 多 种 采样 器 ， 除 了 本 文 提 
到 的 两 种 ， 还 有 CPnest" , emcee 等 。 测 试 它们 对 于 同一 波形 的 后 验 分 布 采样 效果 ， 有 助 
于 模型 的 改进 和 完善 。 我 们 将 在 接 下 来 的 工作 中 进行 深入 研究 。 
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The Detection of Gravitational Waves from Black Hole 


Encounter 


GUO Wei-chang-feng 


(Department of Astronomy, Beijing Normal University, Beijing 100875, China) 


Abstract: Globular clusters and galactic nuclei contain very dense stellar-mass black holes 
(BHs), and when two BHs approach each other, they could form highly eccentric orbits and 
emit gravitational waves (GWs). The detection of such signals will contribute to further 
research on BH physics and enhance our understanding of BH populations. However, such 
a GW source has not yet been detected and has not been widely considered in the curren- 
t study. Therefore, it is important to research its detectability. Based on the numerical 
relativistic simulation, a recent study has given accurate waveforms for parabolic BH en- 
counters. Using the GW burst search tool Minke, a large-scale data set is produced in the 
two-parameter spaces of total mass and distance, and the detection range was quantitatively 
calculated. The result shows that BH encounter is a potential candidate for current ground- 
based detectors. Finally, a focus question that how to distinguish the BH capture from a 
binary BH merger is discussed. Using the posterior distributions of the signal under the two 
samplers, a distinguishing method is proposed. By comparing different posteriors, it is found 
that there is a significant difference in the merger time to of the BH capture and binary BH 


merger signal, which contributes to distinguishing the two GW sources. 


Key words: gravitational wave; stellar-mass black hole; binary black hole 


